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Abstract 

The interpretation of the charged (CCQE) and neutral (NCE) current quasi 
elastic events collected by the MiniBooNE collaboration involves a number 
of unresolved issues. While it has been suggested that the data can be 
explained in terms of an effective nucleon axial mass, M4, the results of our 
theoretical calculations suggest that the CCQE and NCE samples cannot be 
described by the same value of Ma- We argue that the disagreement between 
theory and data may arise from the uncertainties associated with the flux 
average procedure. We also analyze the role of the strange quark in NCE 
interactions and find that, due to a cancellation between proton and neutron 
contributions, it turns out to be negligible. 
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1. Introduction 

The MiniBooNE collaboration has recently collected an extensive data 
set of quasielastic neutrino nucleus scattering events, in both the charged- 
current (CCQE) Q and neutral current (NCE) Q channels, using a Carbon 
target. 

In the CCQE channel, quasielastic neutrino-nucleon interactions are de- 
scribed in terms of the vector form factors Ff ,n (Q 2 ) and F^iQ 2 ) (Q 2 = —q 2 , 
q being the four-momentum transfer, while the superscripts p and n corre- 
spond to proton and neutron, respectively), that have been precisely mea- 
sured in electron-proton and electron-deuteron scattering experiments j3[, 
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and the axial form factor Fa(Q 2 ) 0, 0, @|. In addition, NCE interactions are 
also affected by the form factors F£ , _F| and F% , arising from strange quark 
contributions 0, SIS) 1Q|- The results of recent experiments 0] indicate that 



F*, F£ are vanishing, whereas the axial form factors Fa and F A are assumed 
to be of dipole form, and their Q 2 -dependence is parametrized in terms of 
the axial mass Ma- 

The measured cross sections turn out to be consistently larger than the 
predictions of Monte Carlo simulations carried out using the relativistic Fermi 
gas (RFG) model of the nucleus and the value of the axial mass resulting 
from the world average of the deuterium data, Ma = 1-03 GeV 0. In order 
to bring the predictions of the RFG model into agreement with the data, 
the authors of Refs. [H, 0] use a significantly larger value of the axial mass, 
Ma ~ 1-35 GeV, and introduce the additional parameter k, meant to improve 
the treatment of Pauli blocking. The K2K collaboration also reported a 
large value of the axial mass, Ma ~ 1.2 GeV, resulting form the analysis 



of its sample of CCQE events Moreover, the best fit to the neutral 

current data is obtained using a non vanishing strange quark contribution 
As, determining the value of F^ at Q 2 — Q. 

It has been suggested that the large value of Ma may be regarded as an 
effective axial mass, modified by nuclear effects not taken into account in the 
RFG model However, the results obtained using more advanced models 
appear to rule out this explanation. In fact, numerical calculations carried 
out using realistic nuclear spectral functions, extensively employed in the 
analysis of electron- nucleus scattering data [l2[ , indicate that reproducin g th e 



CCQE measured cross sections requires an even larger value of Ma \L2i, uA 



The purpos e of this work is the extension of the spectral function approach 



of Refs. |13l . Il4j to the description of NCE interactions and the quantitative 
analysis of the Ma and As dependence of the resulting cross sections. The 
main elements of our approach are outlined in Section |2j while the numerical 
results are discussed in Section [3j Finally, in Section [4] we summarize our 
findings and state the conclusions. 

2. Formalism 

We consider the neutral current process 

v,+ 12 C -*v, + X , (1) 



2 



in which a neutrino carrying initial four-momentum k = (E u , k) scatters 
off a Carbon target to a state of four-momentum k' = (E' u ,\s!), the target 
final state being undetected. In the impulse approximation (IA) scheme 



15l |. stating that when the magnitude of the momentum transfer |q| is large 
enough (i) the target nucleus is seen by the probe as a collection of individual 
nucleons and (ii) in the final state the knocked out nucleon and the recoiling 
nucleus evolve independently of one another, the differential cross section can 
be written in the form 

da IA = J d 3 p dE P(p, E) da e i em , (2) 

where da e i em is the neutrino-nucleon cross section and P(p, E) is the spectral 
function of the target nucleus, yielding the probability distribution of finding 
a nucleon of momentum p and removal energy E in the nuclear target. 

2.1. Neutrino nucleon cross section 

The NCE neutrino nucleon cross section in the center of mass frame reads 

da elem _ \M\ 2 ( E'\ ^ 



dVt 64tt 2 (^ + E p f \E V 

where E p is the nucleon energy and is Feynman's invariant amplitude, 
averaged over the spins of the initial state particles and summed over the 
spins of the particles in the final state. 
Feynman's amplitude can be written as 

M = *G F 9^)^(1 —y B )u(k) < N(p')\r z \N(p) >, (4) 
2a/2 v v ' v v ' 

leptonic current hadronic current 

where and v(k') are the Dirac spinors associated with the initial and 
final state neutrino, respectively, the kets \N(p) > and \N(p') > represent the 
initial and final nucleon state, and Jz is the hadronic neutral current. While 
the leptonic current has a simple V-A structure, completely determined by 
the leptons kinematics, the hadronic current is more complex, on account of 
the strong interactions occurring between the nucleon constituents. 

The hadronic neutral weak current can be written in the general form 

<N(p')\r z \N{p) >= 

<N'\ 7 ^(Q 2 ) + ^^F|(g 2 )+ 7 V^l(Q 2 ) \n>, 
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where F±(Q 2 ), F^iQ 2 ) and F^(Q 2 ) are the Dirac, Pauli and axial form factors 
for neutral current interactions, respectively, taking into account the strange 
quark content of the nucleon. 

The form of the neutral weak current 

Jz = J-z- 2sin 2 8 w J em , (6) 

where J3, J em and 9w are the third component of the isospin current, the 
electromagnetic current and Weinberg's angle, respectively, suggests the fol- 
lowing parametrization of the form factors 



F^{Q 2 ) 


1 
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(7) 



where 

F l {Q 2 ) = F?{Q 2 )-F?{Q 2 ) ,' = 1,2, (8) 

and s indicates the strange quark contribution. As stated above, F( and _F| 
are vanishing, while F A is assumed to have a dipole Q 2 dependence 

n(Q 2 ) = . A ' 5 , (9) 

As being the strange quark contribution to the nucleon spin at Q 2 = 0. 

Following Ref. [8J, we parametrize Feynman's amplitude Ai in terms of 
six contributions according to 

\M\ 2 = 4G 2 F (V U + V 12 + V 22 + A + V M + V A2 ) , (10) 
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with 



V n = 4(Ff) 2 [p ■ kk' -p' + p' -kk' -p- M 2 k ■ k'] , 
V 12 = —AF[F% k ■ k'(p' -p)-(k- k'), 



22 



2(^|) 2 
M 2 



k-k' [p-kp' -k+p- k'p' ■ k' + M 2 k ■ k'] 



AD 



A = 4{G A ) 2 [p ■ kk' -p' +p' -kk' -p + M 2 k ■ k'] 
V A i = 8G A F* [p-kp' -k' -k- p'p ■ k'] , 
V A2 = AGaF* k ■ k'(k + k') ■ (p + p'), 



where M is the nucleoli mass and Ga 
9a = F A (Q 2 = 0). 



-FX(g A ~ Ml {9 a + As), with 



2.2. Target spectral function 

Accurate ab initio calculations of the spectral function P(p, E), based on 
realistic nuclear hamiltonians, can only be carried out for the lightest nuclei 
(A < 4) H Il7l. Il8l . Il9l l20l . |21| and in the limit of uniform nuclear matter 
(A 00) 
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In the case of medium-heavy nuclei, the calculation of 
P(p, E) involves severe difficulties, and one has to resort to some simplifying 
assumptions. 



Within the RFG model [24J, |25| the nucleus is described as a degener- 
ate gas of non-interacting nucleons. According to this picture the spectral 
function takes the simple form 

/Qtt 2 A\ 

PRFGuip, E) = y—j ©(Pf - P) 5(E p -E B + E), (12) 

where E p = a/ M 2 + | p | 2 is the energy of a free nucleon carrying momentum 
p. The Fermi momentum pp and the average binding energy Eb are the 
model parameters, to be adjusted to reproduce the experimental data. 

The spectral function in Eq. [12] is non vanishing only at |p| < pp. 
However, electron-nucleus scattering experiments have provided unambigu- 
ous evidence of strong nucleon-nucleon correlations, that give rise to virtual 
scattering processes leading to the excitation of nucleons to states of large 
momentum and removal energy 12|. Hence, the quantitative analysis of 



neutrino-nucleus interactions requires a more realistic spectral function, tak- 
ing into account correlation effects. 
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In our work we have used the Carbon spectral function of Ref. [261], ob- 
tained within the Local Density Approximation (LDA) combining the infor- 
mation extracted from measurements of the coincidence (e, e'p) cross section 
with theoretical calculations of the spectral function of uniform nuclear mat- 
ter at different densities. 

The resulting P(p, E) consists of two contributions (26| 

Plda{p, E) = P M f(p, E) + P corr (p, E) , (13) 

arising from the nuclear mean field and from nucleon-nucleon correlation. 
The mean field spectral function reads 

Pmf(p,E) = J2 Z n\Mp)\ 2 F n (E-E n ) , (14) 

n 

In the above equation, n (p) is the squared momentum-space wave function 
of the n-th shell model state, whose width is described by the Lorentzian 
F n (E—E n ), Z n is the corresponding spectroscopic factor and the sum extends 
to all states belonging to the Fermi sea. In the absence of correlation Z n — > 1 
andP„(P-P„) -+ 6(E - E n ). 

The correlation contribution to the LDA spectral function can be written 
in the form 

P corr ,(p, E) = J d 3 r p A (r) P^(p, E- p = p A (r)) , (15) 

where Pa(?) is the nuclear density profile and P c ^ r A /(p, E; p) is the correlation 
part of the nuclear matter spectral function at density p, whose calculation 



is described in Ref. [26(. Correlation effects turn out to be sizable, leading 



20% of the strength to the region of large momentum (|p| > pp) and large 



26 



energy 

In the IA scheme statistical correlations leading to the suppression of the 
phase-space available to the final state nucleon, generally referred to as Pauli 
Blocking (PB), are not taken into account. In order to introduce their effect 
in our calculations, we have modified the spectral function according to [3] 

P(p, E) => P(p, E) 0(|p + q| - p F ) , (16) 

where q is the momentum transfer and pp is the average Fermi momentum 
of the nucleus, defined as 

p F = d 3 rp A (r)p F (r) , (17) 
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Figure 1: NCE flux averaged Q 2 -distribution f or different values of the axial mass. The 
data points are taken from Ref. 0] 



with 

( ^ 2 p A {r ^ 1/3 
PF{r) = I — 

For a Carbon target, Eqs. (ITT)) and (fTHI) lead to = 225 MeV. The inclusion 
of PB, while leaving unaffected the cross sections at large Q 2 , leads to an 
appreciable quenching in the region of low Q 2 [loT ]. 

3. Results 

We have computed the Q 2 -distribution, avera ged over the MiniBooNE 



flux, using the Carbon spectral function of Ref. 26| . Figure [Tj shows the 
results corresponding to different values of the axial mass and As = 0, com- 
pared to the experimental data of Ref. j^] . It clearly appears that the value 
of the axial mass yielding a good fit of the MiniBooNE CCQE distribution, 



Ma = 1.6 GeV [14J, does not reproduce NCE data. 

To illustrate the dependence of our results on As, in Fig. [2] we show the 
flux averaged Q 2 -distribution for neutrinos interacting with a Carbon target, 
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Figure 2: NCE flux averaged (^-distribution for different values of As. 

obtained using Ma = 1.03 GeV and setting As = and As = —0.19, the 
latter being the lowest value that can be found in the literature BOSH- 



It is apparent that the distribution is nearly independent of As. As a 
consequence, the results displayed in Fig. CD have been obtained neglecting 
the strange quark contribution to the axial form factor. 

In order to analyze the difference between the proton and neutron contri- 
butions, in Fig. [3] we show the same distributions as in Fig. [2j calculated for 
a A = 12 target consisting of protons or neutrons only. It turns out that the 
strange quark produces a suppression of the neutron contributions and an 
enhancement of the proton contribution of about the same size. As a result, 
the two effects largely cancel each other in the Carbon Q 2 -distribution, as 
seen in Fig. [2j 

The role of strange quarks had been previously discussed in Refs. (8l lijj]. 
whose authors proposed to determine As from the ratio 

, x NCE / , \ NCE 

w) <{w) • (19) 

^ / neutron \ ^ / proton 

that does not suffer from the uncertainties associated with the incoming 
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Figure 3: NCE proton and neutron contributions to the Carbon Q 2 -distribution. 



neutrino flux and is very sensitive to variations of As. 



4. Conclusions 

Our work indicates that the theoretical analysis of the MiniBooNE NCE 
data sample involves the same difficulties already emerged in the studies of 
CCQE interactions 
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The results discussed in Section [3j showing that it is impossible to de- 
scribe both the CCQE and NCE data sets using the same value of the axial 
mass, confirm that nuclear effects not included in the oversimplified RFG 
model cannot be taken into account through a modification of Ma- In this 
context, it has to be pointed out that the need of a larger Ma to reproduce 
the measured NCE Q 2 -distribution is not likely to be ascribable to different 
nuclear effects in the CCQE and NCE channels. In fact, the ratio between 
the Q 2 -distributions obtained from the RFG model and the spectral function 
approach, providing a measure of the effects of nuclear dynamics, turns out 
to be nearly identical for CCQE and NCE. The difference does not exceed 
2% over the whole Q 2 range. 
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Our analysis also shows that the strange quark contribution to the cross 
section of nuclei with equal number of protons and neutrons is vanishingly 
small. As a consequence, the possibility of improving the agreement between 
MC simulations and Carbon data adjusting the value of As appears to be 
ruled out. 



The authors of Ref . [14J argued that the disagreement between theory and 
MiniBooNE CCQE data may be due to the uncertainties associated with the 
flux average procedure, as the resulting cross section at fixed energy and 
scattering angle of the outgoing muon picks up contributions from different 
kinematical regions, where different reaction mechanisms are known to be 
dominant. 

This uncertainty also affects the flux averaged NCE differential cross sec- 
tion, which is given in bins of reconstructed Q 2 0], defined as 

Q 2 rec = 2MT = 2M^2 T i , (20) 

i 

where M is the nucleon mass and T is the sum of the kinetic energies of the 
final state nucleons. 

In order to provide results that can be compared to data in a meaningful 
fashion, theoretical models of neutrino-nucleus interactions must be based on 
a consistent description of the broad kinematical range corresponding to the 
relevant neutrino energies. In the quasi elastic sector, this amounts to taking 
into account, besides single- nucleon knock out, multi- nucleon knock out as 
well as processes involving the nuclear two-body currents, whose contribution 



is expected to be significant [29 
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